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Abstract 

The phase diagram of the ferromagnetic Kondo model for manganites is investigated using computational tech- 
niques. In clusters of dimensions 1 and 2, Monte Carlo simulations in the limit where the localized spins are classical 
show a rich low temperature phase diagram with three dominant regions: (i) a ferromagnetic phase, (ii) phase separa- 
tion between hole-poor antiferromagnetic and hole-rich ferromagnetic domains, and (iii) a phase with incommensurate 
spin correlations. Possible experimental consequences of the regime of phase separation are discussed. Studies using 
the Lanczos algorithm and the Density Matrix Renormalization Group method applied to chains with localized spin 
1/2 (with and without Coulombic repulsion for the mobile electrons) and spin 3/2 degrees of freedom give results in 
excellent agreement with those in the spin localized classical limit. The Dynamical Mean Field (D — oo) approxima- 
tion was also applied to the same model. At large Hund coupling phase separation and ferromagnetism were identified, 
again in good agreement with results in low dimensions. In addition, a Monte Carlo study of spin correlations allowed 
us to estimate the critical temperature for ferromagnetism T^ M in 3 dimensional clusters. It is concluded that T^ M 
is compatible with current experimental results. 



1 Introduction and Main Results the development of La-manganite sensors is still at a very- 
early stage, a more fundamental approach to the study 
Materials that present the phenomenon of "colossal" f manganites is appropriate and, thus, theoretical guid- 
magnetoresistance are currently under much experimen- ance i s needed. The existence of correlation effects in 
tal investigation due to their potential technological ap- the fairly dramatic magnetic, transport, and magneto- 
plications. Typical compounds that have this phe- transport properties of doped La-manganites reinforces 
nomenon are ferromagnetic (FM) metallic oxides of this notion, 
the form Ri_ x X x Mn0 3 (where R = La,Pr,Nd; X = Sr, 

Ca, Ba, Pb) jjl Sj. As an example, a decrease in resistiv- The early theoretical studies of models for mangan- 
ity of several orders of magnitude has been reported in ites concentrated their efforts on the existence of fer- 
thin films of Ndo.7Sro.3Mn03 at magnetic fields of 8 Tes- romagnetism. The so-called "Double Exchange" (DE) 
las H . The relative changes in resistance for the mangan- model ^| explained how carriers improve their kinetic 
ites can be as large as AR/R ~ 100, 000%, while in mag- energy by forcing the localized spins to become ferromag- 
netic superlattices Co/Cu/Co the enhancement is about netically ordered (this phenomenon is quite reminiscent 
100%. This result suggests that manganites indeed have of the Nagaoka phase discussed in models for cuprates). 
technological potential since large changes in resistance However, in spite of this successful explanation of the ex- 
can be obtained at fixed temperature upon the applica- istence of ferromagnetism at low temperature several fea- 
tion of magnetic fields, opening an alternative route for tures of the experimental phase diagram of manganites 
next generation magnetic storage devices. However, since remain unclear, and they are likely beyond the DE model. 



1 



Actually, the phase diagram of Lai_ x Ca x Mn03 is very 
rich with not only ferromagnetic phases, but also regions 
with charge-ordering and antiferromagnetic correlations at 
x > 0.5 ||, and a poorly understood "normal" state above 
the critical temperature for ferromagnetism, T C FM , which 
has insulating characteristics at x ~ 0.33. Finding an in- 
sulator above T™ is a surprising result since it would 
have been more natural to have a standard metallic phase 
in that regime which could smoothly become a ferromag- 
netic metal as the temperature is reduced. Some theories 
for manganites propose that the insulating regime above 
Tf M is caused by a strong correlation between electronic 
and phononic degrees of freedom . Other proposals in- 
clude the presence of Berry phases in the DE model that 
may lead to electronic localization || . On the other hand, 
the regime of charge ordering has received little theoretical 
attention and its features remain mostly unexplored. To 
complicate matters further, recent experiments testing the 
dynamical response of manganites have reported anoma- 
lous results in the ferromagnetic phase using neutron scat- 
tering [^| , while in photoemission experiments ]Ic| the pos- 
sible existence of a pseudogap above the critical tempera- 
ture was reported. 

The rich phase diagram of the manganites described 
above, plus the several experimental indications of strong 
correlations in the system, deserves a systematic theoreti- 
cal study using state-of-the-art computational techniques. 
These methods are unbiased and can provide useful infor- 
mation on models for manganites in a regime of couplings 
that cannot be handled perturbatively or exactly. How- 
ever, the large number of degrees of freedom and associ- 
ated couplings of a full Hamiltonian model for mangan- 
ites makes this approach quite cumbersome. In principle 
the two e g active orbitals per Mn-ion must be included, 
in addition to the t% g electrons. Also phonons should be 
incorporated to fully describe these materials. However, 
as a first step towards a theoretical understanding of the 
behavior of models for manganites, in this paper it was 
decided to work only in the electronic sector (i.e. leaving 
aside phonons) and with just one orbital per site (i.e. keep- 
ing only one e g orbital of the two available, which in prac- 
tice amounts to assuming a static Jahn- Teller distortion) . 
In addition, the ti g degrees of freedom are here assumed 
to be localized i.e. no mobility is given to these electrons. 
They basically provide a spin background in which the 
e g electrons move, with a Hund term that couples all ac- 
tive electrons per Mn-ion. Under these assumptions a rich 
phase diagram was obtained, as explained in detail in the 
rest of the paper. It is left for future work the analysis of 
the influence of phonons and orbital degeneracy into the 
fairly complicated phase diagrams reported here. 

Some of the main results of the present effort are sum- 
marized in Figs.l and 2 where the phase diagrams of the 
ferromagnetic Kondo model, using classical spins to repre- 
sent the ti g spins, are presented for low dimensional sys- 



tems. These figures were reported previously in a short 
version of this paper pd| | but here considerable more de- 
tails, as well as a large set of novel results, are provided. 
Three dominant regimes have been identified: (1) a ferro- 
magnetic region in excellent agreement with the DE mech- 
anism, (2) phase separation between hole-rich ferromag- 
netic and hole-poor antiferromagnetic regions, and (3) an 
intermediate phase with short-range incommensurate cor- 
relations. The regime of phase separation was previously 
conjectured to exist in this model from studies of Hamilto- 
nians for manganites at large Hund coupling in one dimen- 
sion fl2]| . It is important to remark that phase separation 
is currently widely discussed in the context of high tem- 
perature superconductors since models for the cuprates, 
such as the t — J model in two dimensions, present densi- 
ties that cannot be uniformly stabilized in a given volume 
by suitably selecting a chemical potential |l3|, [l4[ [l5| . Af- 
ter the introduction of 1 jr Coulombic interactions the re- 
sulting hole-rich regions become unstable against the for- 
mation of "stripe" configurations |l6|, [Tt]] . If the tenden- 
cies towards phase separation reported in this paper are 
realized in the manganites, a similar phenomenon may 
likely occur i.e. neutron scattering experiments could re- 
veal evidence of stripe configurations in compounds such 
as Lai_ x Ca x Mn03 as it occurs in the cuprates Jlq ]. The 
pseudogap features found in photoemission JTo| could also 
be related to this phenomenon. 
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Figure 1: Phase diagram of the ferromagnetic Kondo 
model in one dimension obtained with Monte Carlo tech- 
niques in the limit where the localized spins are classical. 
PS, FM, and IC denote regions with phase separation, fer- 
romagnetism, and incommensurate correlations. For de- 
tails see the text. 
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Figure 2: Same as Fig.l but for the case of two dimen- 
sions. The boundary between the PS and IC regimes was 
difficult to obtain numerically and, thus, it is not sharply 
defined in the figure (gray region). 



The bulk of the paper is devoted to the discussion of 
the numerical evidence that provides support to these pro- 
posed phase diagrams. Results were not only obtained in 
dimensions 1 (Sees. Ill and IV) and 2 (Sec.V), but also 
in 3 (Sec. VII) and oo (Sec. VI) using a variety of numer- 
ical techniques. Both classical and quantum mechanical 
iig degrees of freedom were analyzed on chains. Exper- 
imental consequences of our results are discussed, spe- 
cially those related with the existence of phase separation 
(Sec. VIII). There is a clear underlying qualitative univer- 
sality between results obtained on lattices with different 
coordination number and using different algorithms. This 
universality lead us to believe that the results reported 
here contain at least the dominant main features of the 
phase diagram corresponding to realistic electronic mod- 
els for manganites that have been widely discussed before 
in the literature, but which were not systematically stud- 
ied using computational methods. 



where c- la are destruction operators for electrons at site i 
with spin a, and Si is the total spin of the ti g electrons, 
assumed localized. The first term is the electron transfer 
between nearest-neighbor Mn-ions, Jh > is the ferro- 
magnetic Hund coupling, the number of sites is L, and the 
rest of the notation is standard. The boundary conditions 
used in the present study are important for some results, 
and they will be discussed later in the text. The electronic 
density of e g electrons, denoted by (n), is adjusted using a 
chemical potential fi. In several of the results shown below 
the spin Si will be considered classical (with |Si| = I). In 
ID both classical and quantum mechanical t2 g spins will 
be studied, the latter having a realistic spin 3/2 but also 
considering spin 1/2 for comparison. Phenomenologically 
Jh ^> t, but here Jn/t was considered an arbitrary pa- 
rameter, i.e. both large and small values for Jn/t were 
studied. Below some calculations were also carried out 
including a large on-site Coulombic repulsion among the 
mobile electrons. 

For a one-dimensional chain (open ends), or for a 
one-dimensional ring with L even and periodic or an- 
tiperiodic boundary conditions (PBC and APBC, respec- 
tively) the model is particle-hole symmetric with respect 
to (n) = 1 by simply transforming c\+ — > (— l) 1 ^ and 

c\, — * — (— l) 1 ^! for the mobile electrons. In this case the 
density is transformed as (n) — > 2 — (n) . Similar transfor- 
mations can be deduced for clusters of dimension larger 
than one. Then, here it is sufficient to study densities 
(n> < 1. 

The FM Kondo model Eq.(l) with classical spins can 
be substantially simplified if the limit Jh — > oo is also 
considered. In this situation at every site only the spin 
component of the mobile electrons which is parallel to the 
classical spin is a relevant degree of freedom. The best way 
to make this reduction in the Hilbert space is by rotating 
the Cio- operators into new operators d- la using a 2 x 2 
rotation matrix such that the transformed spinors point 
in the direction of the classical spin. Explicitly, the actual 
transformation is: 

c iT = cos(0i/2) d iT - sm{9 i /2)e- l4 ' i d n , (2.1) 

c u = sm(0i/2)e i?ii d iT + cos(6» ; /2) d n . (2.2) 

The angles 9i and (pi define the direction of the classical 
spin at site i. After this transformation the new Hamilto- 
nian becomes 



2 Model, Symmetries and Algo- 
rithm 



Hj H=00 = - } (tm.nd^dn^ + h.C.), (3) 
<m,n) 



where the down component of the new operators has been 
The ferromagnetic Kondo Hamiltonian H [L9| studied in discarded since the Jh — > oo limit is considered. The 
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i.e. it is complex and dependent on the direction of the 
classical spins at sites m,n |y. In the limit Jh = oo the 
band for the transformed spinors with up spin has itself 
a particle-hole symmetry, which exists for any arbitrary 
configuration of classical spins. This can be shown by 
transforming — ► (— l) 1 ^, and noticing that after this 
transformation the resulting Hamiltonian matrix is simply 
the conjugate of the original. Then, the eigenvalues are 
unchanged, but (dLdj-f) — > 1 — (dLdif). 

The partition function of the FM Kondo model with 
classical spins can be written as 

Z = n(J Q <*0iSin0ijT d<f>^J tr c (e-^). (5) 

For a fixed configuration of angles {9i,cj>i} the Hamilto- 
nian amounts to non-interacting electrons moving in an 
external field. This problem can be diagonalized exactly 
since the Hamiltonian is quadratic in the fcrmionic vari- 
ables. This diagonalization is performed simply by calling 
a library subroutine in a computer program. If the 2L 
eigenvalues are denoted by e\ the resulting partition func- 
tion becomes 

Z = Il{j a ddiS ' m0i J o #ij II( 1 + e " /3eA )- ( 6 ) 

The integrand in Eq.(6) is positive and, thus, a Monte 
Carlo simulation can be performed on the classical spin an- 
gles without "sign" problems. This was the procedure fol- 
lowed in our study below. The number of sweeps through 
the lattice needed to obtain good statistics varied widely 
depending on the temperature, densities, and couplings. 
In some cases, such as in the vicinity of phase separation 
regimes, up to 10 6 sweeps were needed to collect good 
statistics. Measurements of equal-time spin and charge 
correlations for the mobile electrons were performed by 
transforming the operators involved using the basis that 
diagonalizes the Hamiltonian for a fixed configuration of 
classical spins. Dynamical studies of the optical conduc- 
tivity <j{u>) and the one particle spectral function A(p, ui) 
could also be performed following a similar procedure, but 
their detailed study is postponed for a future publication. 
The analysis of the FM Kondo model in the case of quan- 
tum mechanical ti g degrees of freedom was performed at 
zero temperature using standard Lanczos |l4], ^(J and Den- 
sity Matrix Renormalization Group (DMRG) plj tech- 
niques. The study at infinite dimension was carried out 
with the Dynamical Mean-Field approach |l!| . 

3 Results in D=l with Classical 
Localized Spins 

In this section the computational results that led us to pro- 
pose Fig.l as the phase diagram of the FM Kondo model in 



one dimension using classical spins are presented. Results 
for quantum mechanical ti g spins are also provided. 

3.1 Ferromagnetism 

The boundary of the ferromagnetic region in Fig.l was 
obtained by studying the spin-spin correlations (between 
the classical spins) defined in momentum space as S(q) = 
(1/L) J2j m e t ( J ' -m ^' 9 (Sj • S m ), using a standard notation, 
at the particular value of zero momentum. L is the num- 
ber of sites. The analysis was performed for couplings 
J H /t = 1,2,3,4,8,12 and 18. For the last four (strong) 
couplings and in the region of densities that were found 
to be stable (see next section for a discussion) the real 
space spin-spin correlations (Sj ■ S m ) have a robust tail 
that extends to the largest distances available on the clus- 
ters studied here. This is obviously correlated with the 
presence of a large peak in S(q) at zero momentum. Fig. 3 
shows representative results at Jn/t = 12 using open 
boundary conditions (OBC) and L = 24. A robust fer- 
romagnetic correlation is clearly observed even at large 
distances. The strength of the tail decreases as (n) in- 
creases in the region (n) > 0.5, and it tends to vanish at a 
density ^0.78 which corresponds to the boundary of the 
ferromagnetic region in Fig.l. A similar behavior is ob- 
served for smaller values of the coupling. It is interesting 
to note that the maximum strength of the spin-spin corre- 
lation tail appears at (n) m 0.50. This result is compatible 
with the behavior expected at Jh — oo where large and 
small densities are exactly related by symmetry (as dis- 
cussed in Sec. II) and, thus, the ferromagnetic correlations 
should peak at exactly (n) = 0.50. Working at Jn/t < 8 
this qualitative behavior is washed out and in this regime 
the correlations at densities (n) < 0.50 are very similar. 

Special care must be taken with the boundary con- 
ditions (BC). Closed shell BC or open BC are needed 
to stabilize a ferromagnetic spin arrangement. If other 
BC are used the spin correlations at short distances are 
still strongly FM (if working at couplings where ferromag- 
netism is dominant), but not at large distances where they 
become negative since a kink appears separating two re- 
gions of opposite total spin, with each region having all 
spins aligned in the same direction. This well-known ef- 
fect was observed before in a similar context |l2], ^2| and 
it does not present a problem in the analysis shown below. 
Actually results using other boundary conditions and lat- 
tice lengths up to 60 sites are compatible with the data 
of Fig. 3. In particular it can be shown that the ferromag- 
netic correlations persist when the lattice size is increased. 
Fig. 4a shows S(q = 0) versus temperature for several lat- 
tice sizes. At small temperature the sum over sites of the 
correlations grow with L showing that the FM correlation 
length is larger than the size of the chains used here. How- 
ever, in ID the Mcrmin- Wagner theorem forbids the exis- 
tence of a finite critical temperature and, thus, eventually 
S(q = 0) must tend to saturate at any fixed finite temper- 
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ature T as L increases, as it already occurs for T > 0.02 
in Fig. 4a. In spite of this subtle detail, using closed-shell 
boundary conditions the numerical data presented in this 
subsection supports the presence of strong ferromagnetic 
correlations in the one dimensional FM Kondo model, and 
it is reasonable to expect that at zero temperature the 
model in the bulk will develop a finite magnetization in 
the ground state. Fig. 4b contains the spin-spin correla- 
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Figure 3: Spin-spin correlations among the localized 
spins (assumed classical) obtained with the Monte Carlo 
method working on a ID system, and at the coupling, tem- 
perature, densities, and lattice size shown. Open bound- 
ary conditions were used. Clear strong ferromagnetic cor- 
relations are observed. 

tions vs distance parametric with temperature. They show 
that even at relatively high temperatures, the correlations 
at short distances are clearly ferromagnetic, an effect that 
should influence on transport properties since carriers will 
react mainly to the local environment in which they are 
immersed. Fig. 5 contains information that illustrates the 
formation of ferromagnetic spin polarons at relatively high 
temperatures. Fig. 5a shows the spin correlations at low 
electronic density. At T — i/10 and distance 1, the corre- 
lation is only about 15% of the maximum. However, if this 
correlation is measured in the immediate vicinity of a car- 
rier using (n^Si ■ S i+ i), where n\ = J^a n i<y( l ~ n i-a) and 
n ia is the number operator for e s -electrons at site i and 
with spin projection <r, then the correlation is enhanced to 
40%. Thus, it is clear that FM correlations develop in the 
vicinity of the carrier [ p3| . This spin polaron apparently 
has a size of 3 to 4 lattice spacings in our studies. 
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Figure 4: (a) Temperature dependence of the spin cor- 
relations in the one dimensional FM Kondo model using 
classical spins and the Monte Carlo method. Shown is the 
zero momentum spin structure factor S(q) (for the classi- 
cal spins) at the density, coupling, and lattice sizes indi- 
cated. Closed shell boundary conditions were used (peri- 
odic boundary conditions (PBC) for L = 10 and 30, and 
antiperiodic boundary conditions (APBC) for L = 20 and 
40); (b) Spin-spin correlations at density (n) = 0.7, using 
a 10 site chain, PBC, and Jjj/t = 8.0. The temperatures 
are indicated. 



5 



: 0.5 



0.0 



n 1 — 


— ■— 
























T/t 








•-- 


-• : 1/75 








o-- 


-O : i/50 








■ -- 


-■ : 1/25 




- \\h 




□ -- 


-a : 1/10 






~~». 






'V 








< 












"~a 














T 














- - O- 


• 




^ u 



0.10 




2 3 
I 



Figure 5: (a) Spin-spin correlations at (n) =0.1 using a 
chain of 10 sites with PBC. The Hund coupling is Jn/t = 
8.0; (b) Similar as (a) but measuring the spin correlations 
near a conduction electron using (n?S, • Si + i), where nf is 
defined in the text. 



3.2 Phase Separation 

One of the main purposes of this paper is to report the 
existence of phase separation in the FM Kondo model. To 
study this phenomenon in the grand-canonical ensemble, 
where the Monte Carlo simulations are performed, it is 
convenient to analyze (n) versus /i. If is discontin- 

uous then there are densities that can not be estabilizcd, 
regardless of the value of fi. The results shown in Fig. 6 
obtained at Jn/t = 12 clearly show that indeed phase sep- 
aration occurs in the FM Kondo model. The discontinuity 
is between a density corresponding to the antiferromag- 
netic regime (n) = 1.0 and <~ 0.77 where ferromagnetic 
correlations start developing, as shown before in Fig. 3. In 
the case of the canonical ensemble these results can be 
rephrased as follows: if the system is initially setup with 
a density in the forbidden band, it will spontaneously sep- 
arate into two regions having (i) antiferromagnetic (AF) 
correlations and no holes, and (ii) FM correlations and 
most of the holes. 

The existence of phase separation can also be deduced 
from the actual Monte Carlo runs since they require large 
amounts of CPU time for convergence in the vicinity of 
the critical chemical potential /x c . The reason is that in 
this regime there are two states in strong competition. 
Qualitatively this effect can be visualized analyzing (n) 
as a function of Monte Carlo time. Fig. 7 shows such a 
time evolution when the chemical potential is fine tuned 
to its critical value \i c ~ — 6.69812i at Jn/t = 8. Wild 
fluctuations in (n) are observed with frequent tunneling 
events covering a large range of densities. A change in /x 
as small as 0.001 or even smaller reduces drastically the 
frequency of the tunneling events, and makes the results 
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Figure 6: Electronic density (n) vs the chemical poten- 
tial (i obtained with the Monte Carlo technique applied 
to the one dimensional FM Kondo model with classical 
spins. Coupling, temperature, lattice sizes, and bound- 
ary conditions are indicated. The discontinuity suggests 
that some densities are unstable, signaling the presence of 
phase separation. 
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Figure 7: Monte Carlo time evolution of the density (n) 
at the particular value of fj, where the discontinuity takes 
place working at Jn/t = 8, T = t/75, using L = 20 sites, 
and PBC. Frequent tunneling events are observed showing 
the competition between two states as in a first order phase 
transition. 



more stable although certainly strong fluctuations remain 
in a finite window near fi c . Fig. 8 illustrates this effect 
showing a histogram that counts the number of times that 
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a density in a given window of density is reached in the 
simulation. As [i crosses its critical value the histograms 
change rapidly from having a large peak close to (n) ~ 1 
to a large peak at (n) ~ 0.75. At densities that are stable 
away from the phase separation region these histograms 
present just one robust peak. 

The qualitative behavior exemplified in Fig. 6 was also 
observed at other values of Ju/t- For instance, Fig. 9a con- 
tains resuls for Ju/t = 4 which are very similar to those 
found at a larger coupling. In a weaker coupling regime 
the discontinuity is reduced and now the competition is be- 
tween the AF state and a state with incommensurate cor- 
relations rather than ferromagnetism. As example, Fig. 9b 
contains results for Jn/t = 2. The discontinuity is lo- 
cated near /i c ~ — 1.35i and the two states competing in 
this regime actually have very similar properties. Analyz- 
ing, for instance, just S(q) would have not provided an 
indication of a sharp discontinuity in the density, unlike 
the case of a large Hund coupling where the peak in the 
spin structure factor jumps rapidly from q = n to at the 
critical chemical potential. 
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Figure 8: Hystogram for the total number of electrons 
obtained at Jn/t = 8, T = t/75, using periodic boundary 
conditions on a 20 sites chain. The chemical potentials 
are shown, (a) corresponds to a case where the average 
density is close to (n) — 1, (b) is at the critical chemical 
potential, and at (c) the electronic density is close to 0.77. 
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Figure 9: Electronic density vs n for the FM Kondo 
model with classical spins in one dimension using Monte 
Carlo methods. A clear discontinuity signals the existence 
of phase separation, (a) corresponds to Ju/t — 4, while 
(b) is for Jh /t = 2. Temperatures, chain sizes, and bound- 
ary conditions are indicated. The apparent second discon- 
tinuity located at fi/t ~ —1.55 in (b) is expected to be just 
a rapid crossover. 

3.3 Incommensurate Correlations 

The tendency to develop a spin pattern with incommen- 
surate characteristics can be easily studied in our calcu- 
lations observing the behavior of S(q) as couplings and 
densities are varied. While at large Ju/t the spin struc- 
ture factor is peaked only at the momenta compatible 
with ferro or antiferromagnetic order, a different result 
is obtained as Jn/t is reduced. Fig. 10 shows S(q) at 
Jn/t =1.0 for a variety of densities. The AF peak close to 
(n) = 1 smoothly evolves into a substantially weaker peak 
which moves away from q = n in the range 0.5 < (n) < 1.0. 
The peak position is close to 2kp = (n) , and since the spin- 
spin correlations for the mobile electrons show a similar 
behavior this is compatible with Luttinger liquids predic- 
tions ^4j. In this density range, and for the lattice sizes 
and temperatures used here, the peak in S(q) is broad 
and, thus, it cannot be taken as an indication of long- 
range incommensurate (IC) correlations but rather of the 
presence of IC spin arrangements at short distances. In the 
region of low densities (n) < 0.4, S(q) is now peaked at 
zero momentum which is compatible with the presence of 
robust ferromagnetic spin-spin correlations at the largest 
distances available in the studied clusters. The transition 
from one regime to the other in the intermediate small 
window 0.4 < (n) < 0.5 is very fast and was not studied 
in detail here, but for a second order phase transition it is 
expected to be continuous. In Fig. 11 results for a slightly 
larger coupling Jn/t = 2 are shown. Here the pattern 
is more complicated since apparently the AF peak does 
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not evolve smoothly into the peak at q ~ 7r/2 observed at 
(n) — 0.733, suggesting an interesting interplay between 
spin and charge. 




manganites was predicted theoretically using a Hartree- 
Fock approximation p5| . Our results are compatible with 
these predictions although, once again, it is not clear if 
the IC pattern corresponds to long-range order or simply 
short distance correlations. More work is needed to clarify 
these issues. Nevertheless, the present effort is enough to 
show that the tendency to form IC spin patterns exists in 
the FM Kondo model at small Jji/t- 

3.4 Influence of a Direct Coupling Among 
the Localized Spins 

The results of Sec. III. B show the presence of phase sep- 
aration near half-filling, but not in the opposite extreme 
of low conduction electron density. This is a consequence 
of the absence of a direct coupling among the localized 
spins in Eq.(l). This coupling may be caused by a small 
hybridization between the ti g electrons. If a Heisenberg 
term J'Xwij) ^ ■ coupling the classical spins is added, 
then at (n) = an antiferromagnetic state is recovered 
similarly as at (n) = 1. This term was already considered 



q/n 



Figure 10: S(q) (Fourier transform of the spin-spin cor- 
relations between classical spins) at Jh = t and T = t/75, 
on a chain with 30 sites and PBC. The densities are in- 
dicated. The positions of the peaks indicate the tendency 
to have incommensurate correlations in the FM Kondo 
model. 
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Figure 11: Same as Fig. 10, but working at Jh = 2t and 
T = t/100, with 30 sites and open boundary conditions 
(OBC). 

Note that the presence of IC correlations in models for 
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Figure 12: Phase diagram in the plane J' /t-(n) obtained 
studying the density vs [i calculated using the Monte Carlo 
technique with classical localized spins. The coupling is 
fixed to J jilt = 8 and the temperature is low. Note that 
the regime of low density and J'/t = is somewhat diffi- 
cult to study due to the influence of van Hove singularities 
which produce a rapid change of the density with ^. 

in the study of the strong coupling version of the Kondo 
model . Although a detailed study of the influence of J' 
on the phase diagram Fig.l is postponed for a future pub- 
lication |2(| , here the effects of this new coupling into the 
existence of phase separation are reported. Following the 
same procedure described in Sec. III. B to obtain unstable 
densities, the phase diagram shown in Fig. 12 was found. 
Note that at J' ^ phase separation occurs at large and 
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small electronic densities. In the latter the separation is 
between electron-rich ferromagnetic and electron-poor an- 
tiferromagnetic regions. Then, the experimental search 
for phase separation discussed below in Sec. VIII should be 
carried out at both large and small hole densities. More 
details will be given elsewhere f2(| . 



4 Results in D=l with Quantum 
Localized Spins 

4.1 Quantum Mechanical t 29 S=3/2 spins 

The use of classical spins to represent the t 2g degrees of 
freedom is an approximation which has been used since the 
early days of the study of manganites j|, || . While such an 
approach seems reasonable it would be desirable to have 
some numerical evidence supporting the idea that using 
spin operators of value 3/2 (denoted by Si) the results are 
similar as those obtained with classical spins. Although 
numerical unbiased calculations with spins 3/2 are diffi- 
cult in dimensions larger than 1, at least this issue can be 
addressed numerically with ID chains using the Lanczos 
and Density Matrix Renormalization Group techniques. 
The Hamiltonian is the same as in Eq.(l) but now with 
quantum mechanical degrees of freedom normalized to 1 
(i.e. replacing Si by Si/(3/2)) to simplify the comparison 
of results against those obtained using classical spins. Cal- 
culating the ground state energy in subspaces with a fixed 
total spin in the z-direction, it is possible to study the 
tendency to have a ferromagnetic state in the model using 
the Lanczos technique. The results indicate that there is 
a robust region of fully saturated ferromagnetism, as indi- 
cated in Fig. 13. The actual boundary of this region agrees 
accurately with the results obtained using classical spins 
shown in Fig.l. This reinforces the belief that S=3/2 and 
classical spins produce similar results, at least regarding 
ferromagnetism. Finite size effects are apparently small 
for the chains accessible with the DMRG method, as ex- 
emplified in Fig. 14a where the ground state energy for 
Jii/t = 6.0 using a variety of chain lengths is presented. 

To study other important features of the phase diagram, 
such as phase separation, the compressibility is needed. 
This quantity is defined as 

_ a N 2 e E(N e + 2, L) + E(N e — 2,L) — 2E(N e , L) 

K =— 4 ' 

(7) 

where E(N e ,L) is the ground state energy corresponding 
to a chain with L sites and N e electrons. If k _1 becomes 
negative at some fixed density (remember that here the 
numerical study is in the canonical ensemble), the sys- 
tem is unstable and phase separation occurs. The DMRG 
results for the compressibility in Fig. 14b for several cou- 
plings led us to conclude that the spin 3/2 model also has 
phase separation in the ground state near half-filling, simi- 
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Figure 13: Phase diagram of the FM Kondo model with 
5 = 3/2 localized t2 g spins obtained with the DMRG and 
Lanczos techniques applied to finite chains as indicated. 
The notation is as in Figs.l and 2, i.e. ferromagnetism, 
incommensurability, and phase separation appear both for 
classical and quantum mechanical localized spins. 
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Figure 14: (a) Ground state energy per site vs density at 
Jff/f = 6.0 for the S = 3/2 FM Kondo model in ID using 
DMRG techniques keeping m = 48 states. Results for a 
variety of chain lengths are shown; (b) Inverse compress- 
ibility vs density for the same model and chain lengths as 
in (a), calculated using several couplings. A negative 1/k 
signals an unstable density. 
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larly as in the case of classical spins. Indeed Fig. 14b shows 
that X/k becomes negative for densities in the vicinity of 
(n) ~ 0.85 and larger, for the range of couplings shown. 
Thus, at least qualitatively there are tendencies to phase 
separate in the same region suggested by the simulations 
using classical localized spins. 

However, a difference exists between results obtained 
with classical and quantum mechanical t<i g degrees of free- 
dom: apparently there is a finite window in density be- 
tween the phase separated and FM regimes in the phase 
diagram of Fig. 13. This window could be a finite size ef- 
fect, but the lack of a strong dependence with the chain 
lengths in the results of Fig. 14a led us to believe that it 
may actually exist in the bulk limit. In addition, pre- 
vious studies using the strong Ju/t coupling version of 
the FM Kondo model have also reported an intermediate 
window between phase separation and FM ]T^] , and the 
analysis below for localized spins 1/2 suggests a similar 
result. Studying the spin of the ground state in this inter- 
mediate region here and in Ref.|l^| it has been observed 
that it is finite i.e. apparently partial ferromagnetism ap- 
pears immediately at any finite stable density in the model 
with spins 3/2, at least working at intermediate and large 
Hund couplings |27j . The transition from phase separation 
to FM appears to be smooth in the spin quantum num- 
ber, which is somewhat reminiscent of the results in the 
classical limit where the tail of the spin-spin correlation 
grows with continuity from small to large as the density 
diminishes in the stable region. 
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Figure 15: S(q) for the quantum FM Kondo model using 
localized spins 3/2 on a chain of 12 sites. The technique 
is DMRG, keeping 48 states in the iterations. Densities 
and coupling arc indicated. In the spin correlations used 
to obtain S(q), full spins 3/2 (i.e. not normalized to 1) 
were used. 



Finally, let us analyze whether incommensurate corre- 
lations exist in the spin 3/2 model as it occurs in the 
classical model. Fig. 15 shows S(q) obtained with DMRG 
working on a chain of 12 sites. Although momentum is not 
a good quantum number on a system with OBC, neverthe- 
less using the same definition of the Fourier transform of 
the spin correlation as in the study of the classical system 
with periodic and antiperiodic boundary conditions (PBC 
and APBC, respectively), qualitative information about 
the tendency to form incommensurate structures can be 
gathered. The results of Fig. 15 obtained at very small 
Jn/t indeed show that strong IC correlations exist in the 
ground state of this model, in agreement with the results 
of Figs. 10 and 11 for classical spins. Then, it is concluded 
that the three dominant features of the phase diagram 
Fig.l (PS, FM, and IC) have an analog in the case of the 
spin 3/2 quantum model. This agreement gives support 
to the believe that the results presented below in this pa- 
per for dimensions larger than one using classical spins 
should be qualitatively similar to those corresponding to 
a model with the proper quantum mechanical t2 g degrees 
of freedom. 

4.2 Quantum Mechanical £ 2 g S=l/2 spins 

For completeness, in this paper the special case of local- 
ized spins 1/2 in ID has also been studied. The anal- 
ysis has relevance not only in the context of manganites 
but also for recently synthesized one dimensional materials 
such as Y2_xCa x BaNi05 which have a mobile and a local- 
ized electron per Ni-ion. This compound has been studied 
experimentally |^8) and theoretically p9| , and upon dop- 
ing interesting properties have been observed including a 
metal-insulator transition. As discussed below, the con- 
clusion of this subsection will be that the results for local- 
ized spins 1/2 are qualitatively similar to those obtained 
with classical and spin 3/2 degrees of freedom, i.e. fer- 
romagnetism, incommensurability, and phase separation 
appear clearly in the phase diagram. The Hamiltonian 
used for this study is as defined in Eq.(l) but now with 
Si replaced by a spin 1/2 operator (not normalized to 1). 
The technique used to obtain ground state properties is 
the finite-size DMRG method on chains with up to 40 
sites, a variety of densities, and typical truncation errors 
around 10~ 5 . Some calculations were also performed with 
the Lanczos algorithm on lattices with up to 12 sites. The 
data obtained with the DMRG and Lanczos methods are 
qualitatively similar. 

The main result is contained in the phase diagram 
shown in Fig. 16. Three regimes were identified pof . 
The region labeled FM corresponds to saturated ferro- 
magnetism i.e. the ground state spin is the maximum. 
This regime was found using the Lanczos technique sim- 
ply searching for degeneracies between the lowest energy 
states of subspaces with different total spin in the z- 
direction. For the IC regime, S(q) was calculated using 
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Figure 16: Phase diagram of the FM Kondo model with 
5 = 1/2 localized states obtained with the DMRG and 
Lanczos methods. The length of the chains is indicated. 
The notation FM, PS, and IC is as in previous figures. 
Note the similarity with the results of Fig. 13, up to an 
overall scale. 



the same definition as in subsection III. A but now with 
spin 1/2 operators instead of classical spins of length 1. 
The results shown in Fig. 17 suggest the presence of incom- 
mensurate correlations at small </#/£, at least at short dis- 
tances. Similarly as for the cases of larger localized spins 
(Figs. 10 and 15), the peak in the spin structure factor 
moves smoothly from q = it near (n) = 1 to q — in the 
FM region. The position of the peak is at 2kp. Our study 
also showed that correlated with this behavior the charge 
structure factor N(q) has a cusp at the same position. It 
is important to clarify that in the IC regime of Fig. 17 at 
low temperatures, the ground state has a finite spin but it 
is not fully saturated. Within the accuracy of our study 
the spin varies smoothly as the couplings and density are 
changed reaching its maximum value at the FM boundary. 

The tendency to phase separate in this model can be 
studied calculating (n) vs \i. In the study of this section, 
carried out in the canonical ensemble where {n) is fixed, 
the procedure to obtain \x involves (i) the calculation of 
the ground state energy for a variety of densities, and (ii) 
the addition of — /xiV to the Hamiltonian, where N is the 
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Figure 17: S(q) for the quantum FM Kondo model using 
localized spins 1/2 on a chain of 20 sites. The technique is 
DMRG, keeping up to 60 states in the iterations. Densities 
and coupling are indicated. In the spin correlations used 
to obtain S(q), full spins 1/2 (i.e. not normalized to 1) 
were used. 



total number operator. As /i varies, different density sub- 
spaces become the actual global ground state. If there are 
densities that cannot be stabilized for any value of /x, then 
such a result is compatible with phase separation in the 
model (for more details see Ref. Results are pre- 

sented in Fig.l8.a,b,c for a chain with 20 sites: at small 
Jh /t = 1 all densities are accessible tuning /i, and the cur- 
vature of the ground state energy Eqs vs (n) is positive. 
However, working at Jn/t = 30 the density (n) = 0.9 
becomes unstable, and now Egs vs (n) has less curva- 
ture. Fig.l8d shows again density vs chemical potential 
but now on a larger chain of 40 sites. Here densities be- 
tween 1 and 0.8 are unstable, in agreement with Fig. 18c. 
A similar method to calculate the region of phase separa- 
tion is to evaluate the inverse compressibility since a neg- 
ative value for this quantity indicates phase separation, as 
explained in Sec. IV. A. Using this procedure once again a 
region of < was identified signaling unstable den- 
sities in the system. From the combination of these type 
of analysis performed for several chains, the boundaries of 
phase separation were estimated as shown in Fig. 16. Al- 
though the error bars are not negligible, the presence of 
phase separation is a robust feature of the calculation and 
it is in excellent agreement with the conclusions of previ- 
ous sections. Then, irrespective of the actual value of the 
spin corresponding to the ti g degrees of freedom the phase 
diagram presents universal features, specially robust fer- 
romagnetism, unstable densities, and short-range incom- 
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Figure 18: Results for the FM Kondo model with 5=1/2 
localized spins using the DMRG method, (a) contains (n) 
vs /i for a chain of 20 sites at Jn/t = 1, showing that all 
available densities are accessible; (b) Ground state energy 
per site vs (n) for Jn/t = 1 and 30. The energies for the 
latter were divided by 5 to make both curves of comparable 
magnitude; (c) Same as (a) but for Jn/t = 30. Now 
(n) = 0.90 appears unstable; (d) Results on a 40 sites 
chain at Jn/t — 28. The densities available between 1.00 
and 0.80 (i.e. 0.95, 0.90, and 0.85) are unstable. 



mensurate correlations. Below in Secs.V and VI, it will be 
shown that this universality can actually be extended to 
include higher dimensional clusters. 

4.3 Influence of an On-site Coulomb Re- 
pulsion 

The on-site Coulombic repulsion among electrons in the 
conduction band has been neglected thus far. Although 
it is not expected to produce qualitative changes in the 
physics described before (since a large Jjj /t prevents dou- 
ble occupancy) it would be desirable to have some indica- 
tions of its quantitative influence on the phase diagrams 
discussed in previous sections. Unfortunately, the addition 
of a Hubbard repulsion U complicates substantially the 
many-body numerical studies. The Monte Carlo method 
in the classical limit can only proceed after the Hubbard 
[/-term is decoupled using Hubbard-Stratonovich (HS) 
variables. If such a procedure is followed the simulation 
would run over both angles and HS degrees of freedom, 
and likely a "sign" problem would occur. Such a cumber- 
some approach will not be pursued here. Instead the large 
U/t limit for the case of localized spins 1/2 will be inves- 
tigated using the Lanczos method, which can be applied 
without major complications to the study of Hubbard-like 
systems. The study will be limited to chains due to re- 
strictions on the size of the clusters that can be analyzed 
numerically. In addition, in the strong coupling limit the 
Hubbard model becomes the t — J model, and, thus, the 
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where c is a destruction fermionic operators which includes 
a projector operator avoiding double occupancy. Both Si 
and Si are spin-1/2 operators representing the spin of the 
localized and mobile degrees of freedom, respectively. The 
rest of the notation is standard. 

Once again studying the energy of the ground state for 
subspaces with different total spin projection in the z- 
direction, the existence of ferromagnetism can be studied. 
The boundary of the region with fully saturared ferro- 
magnetism obtained using clusters with 8 and 10 sites is 
shown in Fig. 19a. Strong ferromagnetic correlations ap- 
pear even for small values for Jn/t, suggesting that as U/t 
grows the tendency to favor ferromagnetism increases, as 
noticed also in Ref. |l^] . It is interesting to observe that in 
the region not labeled as FM in Fig. 19a a finite spin exists 
in the ground state (see Fig. 19b), in excellent agreement 
with results for the U/t — case reported in Sec. IV. A and 
IV. B for quantum mechanical t^ g degrees of freedom (see 
also Ref. ||l2f ) . In addition, a tendency towards incom- 
mensurate correlations is observed in this phase, also in 
agreement with previous results (see Fig. 19c). Regarding 
the issue of phase separation, the results of Fig. 16 sug- 
gests that this regime should appear for densities between 
(n) = 1.0 and (n) = 0.90. Unfortunately these densities 
are not accessible on the clusters studied in this subsec- 
tion, and, thus, the confirmation of the existence of phase 
separation at large U ft will still require further numerical 
work ]32|| . Nevertheless, Figs. 16 and 19a, which contain 
results with and without the Coulomb interaction, have 
clear qualitative common trends. Regarding the quanti- 
tative aspects, the Coulomb interaction changes substan- 
tially the scales in the phase diagram specially regarding 
ferromagnetism which now appears at smaller values of 
J H /t @. 

5 Results in D=2 

In this section the computational results that led us to pro- 
pose Fig. 2 as the phase diagram of the FM Kondo model 
in two dimensions are presented. Results in 2D are partic- 
ularly important since manganite compounds with layered 
structure have been recently synthesized |34|. 



5.1 Ferromagnetism 

The search for ferromagnetism in the ground state of the 
2D clusters was carried out similarly as in ID for the case 
of classical t<i g spins. The real space spin-spin correlations 
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(between those classical spins) was monitored, as well as 
its Fourier transform S(q) at zero momentum. Couplings 
Jn/t — 1, 2, 4, 8, and 16 were particularly analyzed. Typ- 
ical results are presented in Fig. 20 where the spin corre- 
lations are presented at a fixed coupling, parametric with 
the electronic density. It is clear that the tails of the corre- 
lations are very robust, and the ferromagnetic correlation 
length exceeds the maximum distance, d maxi available on 
the 6x6 cluster. Plotting the spin correlation at d max 
vs (n), an estimate of the critical density for ferromag- 
netism can be obtained. Combining results from a variety 
of clusters and boundary conditions, the FM boundary of 
Fig. 2 was constructed. Note that for this analysis the use 
of open boundary conditions seem to be the optimal i.e. 
using other boundary conditions kept the ferromagnetic 
character of the system at short distances but modify the 
correlations at large distances, as it occurs for non-closed- 
shell BC in ID. 
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Figure 19: Results corresponding to the FM Kondo model 
with localized spins 1/2 and using the t — J model for 
the mobile electrons (Eq.(8)), obtained with the Lanczos 
method. The boundary conditions were such that a fully 
saturared ferromagnetic states is stable. J/t = 0.2 was 
used, (a) shows the boundary of the FM region using 8 and 
10 sites clusters. IC correlations were observed at small 
JH/t in the non- fully saturated ferromagnetic region; (b) 
the spin of the ground state as a function of Jh for the 
case of 2 holes on the 8 sites cluster; and (c) S(q), the 
Fourier transform of the spin-spin correlations between the 
localized spins, vs momentum for the case of a 10 sites 
cluster. The densities and couplings are indicated. 



Figure 20: Spin-spin correlations (among the classical 
spins) vs. distance obtained at Jn/t = 8, T = t/50, and 
using periodic boundary conditions on a 6 x 6 cluster. Full 
circles, open circles, full squares, open squares, and full 
triangles correspond to densities (n) = 0.807, 0.776, 0.750, 
0.639, and 0.285, respectively. 

The influence of the lattice size can be estimated by 
studying S(q) vs temperature for several clusters. Repre- 
sentative results are shown in Fig. 21. They were obtained 
at Jn/t = 16 and a density close to the PS-FM bound- 
ary. At low temperature S(q = 0) clearly grows with the 
lattice size due to the strong ferromagnetic correlations. 
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Figure 21: S(g) at zero momentum vs. temperature, 
working at (n) = 0.72 ± 0.03, Jujt = 16, using open 
boundary conditions, and for the two-dimensional cluster 
sizes shown. 



However, as in the case of ID, the Mermin- Wagner the- 
orem forbids long-range fcrromagnctism in 2D at finite 
temperature and, thus, S(q) should converge to a finite 
constant if the lattice sizes are further increased beyond 
those currently accessible, working at a fixed finite tem- 
perature. The verification of this subtle detail is beyond 
the scope of this paper, and should not confuse the read- 
ers: the presence of very strong ferromagnetic correlations 
at low temperature in the FM Kondo model is clear in the 
present numerical study and it is likely that small cou- 
plings in the direction perpendicular to the planes would 
stabilize ferromagnetic order at a finite temperature. 
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Figure 22: (n) vs ^ on two-dimensional clusters and 
at temperature T = t/50, to illustrate the presence of 
phase separation in the 2D ferromagnetic Kondo model 
with classical spins, (a) contains results at Jn/t = 8 and 
a variety of clusters and boundary conditions, (b) cor- 
responds to Jn/t — 4, a 6 x 6 cluster, and using both 
periodic and open boundary conditions. 



naling the presence of phase separation in the low tem- 
perature regime of the FM Kondo model. Calculating the 
discontinuity in (n) for several couplings the boundary of 
the phase separated regime in the 2D phase diagram of 
Fig. 2 was established. Note that the scales in computer 
time needed to achieve convergence are very large near the 
critical chemical potential where frequent tunneling events 
between the two minima slow down the simulations. For 
Jn/t < 4, it becomes difficult to distinguish between an 
actual discontinuity in (n) and a very rapid crossover and, 
thus, in the 2D phase diagram (Fig. 2) the boundary of 
phase separation at small Hund coupling is not sharply 
defined. 



5.2 Phase Separation 

The computational analysis carried out in this subsection 
show that the phenomenon of phase separation occurs not 
only in one dimension but also in two dimensions (and 
higher). This result indicates that the unstable densities 
found in Sees. Ill and IV are not a pathology of ID clusters 
but its existence is generic of the FM Kondo model. 

Typical numerical results in 2D clusters at low temper- 
ature are shown in Fig. 22a at Jh = St. Data from a 
variety of cluster sizes and boundary conditions are pre- 
sented. Although the results corresponding to the lowest 
density in the discontinuity in (n) are somewhat scattered 
due to size effects, the existence of such discontinuity is 
clear from the figure. Fig. 22b contains similar results but 
now for Ju/t = 4. Once again, (n) is discontinuous sig- 



5.3 Incommensurate Correlations 

The regime of small coupling Jn/t is not ferromagnetic, 
according to the behavior of S(q) at zero momentum, and 
does not correspond to phase separation since all densi- 
ties are stable. It may occur that robust incommensurate 
spin correlations exist here, as it occurs in ID. S(q) is 
presented in Fig. 23 for two representative couplings and 
a large range of densities. The antiferromagnctic peak at 
(it, tt) is rapidly suppressed as (n) decreases, and the posi- 
tion of the maximum in S(q) moves along the (tt, tt)-(tt, 0) 
line (and also along (tt, tt)—(0, tt) by symmetry). The in- 
tensity of the peak is rapidly reduced as it moves away 
from (tt, tt), and the IC character of the correlations is ap- 
parently short-range. However, the numerical study of IC 
phases are notoriously affected by lattice sizes, and thus at 
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at (a) J H /t = 1 and T = t/25, and (b) J H /t = 2 and 
T = t/50. In both cases open boundary conditions were 
used. The densities are indicated. 



this stage it can only be claimed that a tendency to form 
IC spin patterns has been detected in 2D clusters, with- 
out a firm statement regarding their short- vs. long-range 
character. As (n) is further reduced, the peak position 
reaches (tt,0) and (0, 7r) with substantial intensity. Re- 
ducing further the density, a rapid change into the ferro- 
magnetic phase was observed. At a larger Ju/t coupling 
such as 3, a more complicated IC pattern was detected 
with regions where S(q) peaked simultaneously at (tt, 0)- 
(0,7r) and (0,0). Further work is needed to clarify the 
fine details of the spin arrangements in this regime, but 
nevertheless the results given here are enough to support 
the claim that a tendency to form incommensurate cor- 
relations exists in the ground state of the 2D FM Kondo 
model l33l. 



6 Results in D=oo 

The existence of phase separation and ferromagnetism in 
the ground state of the FM Kondo model can also be 
studied in the limit of D — 00. The Dynamical Mean 
Field equation jl9| is solved iteratively starting from a 
random spin configuration, and as a function of tempera- 



ture and density three solutions have been observed hav- 
ing AF, FM, and paramagnetic character. Efforts were 
concentrated on a particular large coupling Jh/W = 4.0 
studying the temperature dependence of the results, where 
W is the half-width of the semicircular density of states 
D(e) = (2/71-1^)^1 - {e/W) 2 for the e g electrons. Par- 
tial results are contained in Fig. 24a. The presence of fer- 
romagnetism at finite doping and antiferromagnetism at 
half- filling are quite clear in the calculations. Close to 
half-filling and at low temperature, the density (n) as a 
function of fj, was found to be discontinuous, in excellent 
agreement with the results already reported in D = 1 
and 2. Fig. 24b provides a typical example obtained at 
T/W = 0.0003 (results at T/W = 0.002 were already pre- 
sented in Ref. (ll)). The phase separation in this figure 
is between antiferromagnetic and ferromagnetic regions. 
However, in Fig. 24b note that at a slightly larger temper- 
ature the separation occurs between hole-poor antiferro- 
magnetic and hole-rich paramagnetic regions. 

For completeness, in Fig. 25 the density of states A{u>) 
for the AF and FM phases is shown at Jh/W = 2 and 
T/W = 0.005 (for details of the calculation see Ref.@). 
The critical chemical potential where the AF and FM 
phases coexist is fi c ~ — 1.40VF. A(u>) in Fig. 25 is cal- 
culated at /i = fi c for both phases. In the two cases the 
density of states splits into upper and lower bands due 
to the large Hund coupling. The width of the upper and 
lower bands is wider for the FM phase, which causes a 
narrower gapped region centered at w ~ 0. Let us now 
consider the process of hole doping starting at (n) = 1 and 
decreasing fj,. In the AF phase at a > a c , the chemical 
potential lies in the gap. However, at u < /i c the chemical 
potential is located already inside the lower band of the 
FM phase, since this band is wider than in the AF phase. 
This suggest that before the lightly doped AF phase is 
realized in the system by decreasing fi, the FM phase is 
instead stabilized. Thus, the discontinuous change from 
the AF to FM phases at /x = /j, c also causes a jump in the 
carrier number. This discontinuity occurs only when the 
bandwidth of the AF phase is considerably narrower than 
that of the FM phase. For this reason phase separation 
exists only in the large Jh region. 

Summarizing, at very low temperature there is a re- 
markable agreement between the D = 00 and D = 1, 2 
results which were obtained using substantially different 
numerical techniques. Such an agreement give us confi- 
dence that the phase separation effect discussed here is 
not pathological of low dimensions or induced by approx- 
imate algorithms but it is intrinsic of the physics of the 
FM Kondo model, and likely it exists in dimension D = 3 
as well. 
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Figure 24: (a) Phase diagram in the D — oo limit work- 
ing at J H /W = 4.0. The "PS(AF-PM)" region denotes 
phase separation (PS) between a hole-poor antiferromag- 
netic (AF) region, and a hole-rich paramagnetic (PM) re- 
gion. The rest of the notation is standard; (b) Density (n) 
vs fi/W obtained in the D = oo limit, Jh/W = 4.0, and 
T/W = 0.0003. The discontinuity in the density is clear. 
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Figure 25: Density of states in the D = oo limit corre- 
sponding to the antiferromagnetic (solid line) and ferro- 
magnetic (dotted line) solutions, working at Jh/W — 2.0 
and T/W = 0.005. The chemical potential is tuned to be 
at its critical value /i c ~ — 1.40W. 

7 Results in D=3 

Results in three dimensions for large enough clusters are 
difficult to obtain with the Monte Carlo algorithm used 
in this paper. The reason is that the CPU time needed 
to diagonalize exactly the problem of an electron moving 
in a fixed spin background grows rapidly with the num- 
ber of sites. Nevertheless, studies using 4 3 clusters carried 
out as part of this project have shown clear indications of 
strong ferromagnetic correlations in a region of parameter 



space compatible with those found in one and two dimen- 
sions where FM dominates, and thus it is reasonable to 
assume that there are ferromagnetic phases in the FM 
Kondo model in all dimensions from 1 to oo. Regarding 
phase separation, the studies on 4 3 clusters cannot provide 
conclusive evidence due to the presence of intrinsic gaps in 
(n) vs fj, caused by size effects. But once again by simple 
continuity between D = 1,2 and D = oo, the existence 
of phase separation in D = 3 is strongly suggested by our 
results. 

In spite of the size limitations of studies in three di- 
mensions, it is possible to obtain useful information about 
the actual value of the critical temperature in the limit 
of Jh = oo i.e. working in a region which corresponds 
to a fully saturated ferromagnetic state at zero tempera- 
ture. In this limit the problem is simplified since for the 
mobile electrons only the spin component in the direction 
of the classical spin survives, as discussed in Sec. II where 
the effective model ("complex double exchange") was de- 
scribed. The absence of a spin index reduces substantially 
the CPU time for diagonalization in the numerical algo- 
rithm, and allowed us to study clusters with 6 3 sites using 
the Hamiltonian Eq.(3). Measuring the spin-spin correla- 
tion (among the classical spins) in real space and reducing 
the temperature, it is possible to study at what temper- 
ature, T* , such correlation becomes nonzero at the maxi- 
mum distance available in a 6 3 cluster. For temperatures 
higher than T* the spin correlations can be accommodated 
inside the cluster and, thus, the ferromagnetic correlation 
length, £ FM , is finite. Using this idea, upper bounds on 
the critical temperature can be obtained using T* > Tf . 
In addition, it is reasonable to assume that in a 3D sys- 
tem the growth of £ FM with T is very rapid once it starts, 
and T* itself may actually provide a good estimate of the 
critical temperature in the bulk. 

Figs.26.a-d contain the spin-spin correlations at Jh = 
oo on the 6 3 cluster parametric with temperature. Re- 
sults for four densities are given, and only some repre- 
sentative temperatures are shown. The spin correlations 
on a 4 3 cluster (not shown) are in good agreement with 
those provided in Fig. 26. Based on this information the 
temperature where the ferromagnetic correlations reach 
the boundary with a nonzero value can be obtained with 
reasonably small error bars. The results are shown in 
Fig. 27a. Once again, assuming a rapid increase of £ FM 
as the temperature is reduced, the results of Fig. 27a can 
be considered as a rough estimate of the actual critical 
temperature. As anticipated in Sec. III. A, T C FM is max- 
imized at (n) = 0.50. In Fig. 27b, S(q) at zero momen- 
tum is presented as a function of temperature for the four 
densities used in Fig. 26. A rapid growth is observed at 
particular temperatures which are compatible with those 
obtained in Fig. 27a using the tail of the real space correla- 
tions. The present results are qualitatively similar to those 
obtained using high temperature expansions (HTE) pa, 
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Figure 26: Spin-spin correlations (among the classi- 
cal spins) vs distance for several temperatures (as indi- 
cated) working on a 6 3 cluster and using the Monte Carlo 
technique, (a) contains results for (n) = 0.50, (b) for 
(n) = 0.34, (c) for (n) = 0.14, and (d) for (n) = 0.08. 



although our estimates for the critical temperatures are 
smaller. For instance, at (n) = 0.50 the HTE prediction 
gives Tf M ~ 0.16f, about a factor 1.5 larger than our 
result. 
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Figure 27: (a) Upper bound for the critical temperature 
vs density deduced from spin correlations on 6 3 clusters, 
working at Jh — oo; (b) S(q — 0)/L vs temperature at 
several densities, obtained with the Monte Carlo method 
on a 6 3 cluster at Jn/t = oo. 



To obtain the critical temperature in degrees, an esti- 
mate of t is needed. Results in the experimental litera- 
ture for the e g electrons bandwidth range from BW ~ 
1 ef (§ to BW ~ 4eV @. Assuming a dispersion 
e p = — 2t(cosp x + cospy + cosp z ), the hopping amplitude 
should be t = BW/12 i.e. between 0.08 and 0.33 eV. 
With this result our estimate for the critical temperature 
is roughly between Tf M - 100 K and 400 K, which are 
within the range observed experimentally. Then, in our 
opinion it is possible to obtain realistic values for T C FM 
using purely electronic models, in agreement with other 
calculations |3(| . 
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8 Conclusions 

In this paper the phase diagram of the ferromagnetic 
Kondo model for manganites was investigated. Using 
a wide variety of computational techniques that include 
Monte Carlo simulations, Lanczos and DMRG methods, 
and the Dynamical Mean Field approach, regions in pa- 
rameter space with (i) robust ferromagnetic correlations, 
and (ii) phase separation between hole-poor antiferromag- 
nctic and hole-rich ferromagnetic regions were identified. 
In addition, incommensurate spin correlations were ob- 
served in dimension 1 and 2 at small Jn/t. The criti- 
cal temperature towards ferromagnetism for the case of 
a three dimensional lattice was also estimated, and the 
results are compatible with experiments for manganites. 
The agreement between the results obtained with differ- 
ent computational techniques, working at several spatial 
dimensions, and using both classical and quantum me- 
chanical localized spins in the case of chains lead us to be- 
lieve that the conclusions of this paper are robust and they 
represent the actual physics of the ferromagnetic Kondo 
model. 

The novel regime of phase separation is particularly in- 
teresting, and possible consequences of its existence in 
manganites can be envisioned. Experimentally, phase sep- 
aration can be detected using neutron diffraction tech- 
niques if the two coexisting phases have different lattice 
parameters as it occurs in La 2 Cu04 + a, a Cu-oxide with 
hole-rich and hole-poor regions |3!| . NMR and NQR spec- 
tra, as well as magnetic susceptibility measurements, can 
also be used to detect phase separation [|l], |4lj since a 
splitting of the signal appears when there are two dif- 
ferent environments for the ions. Note also that in the 
regime where AF and FM coexist S(q) presents a two peak 
structure, one located at the AF position and the other at 
zero momentum. This also occurs in a canted ferromag- 
netic state and, thus, care must be taken in the analysis 
of the experimental data. Actually, recent experimental 
results by Kawano et al. [ E| are in qualitative agreement 
with the results of Fig. 24a since these authors observed 
a reentrant structural phase transition accompanied by 
"canted ferromagnetism" below Tf M , at 0.10 < x < 0.17 
in Lai_ x Sr x Mn03. In addition the polaron-ordered phase 
reported by Yamada et al. |E§ can be reanalyzed using 
the results of this paper since it is known that the AF 
phase in 3D manganites is orthorhombic while the FM 
is pseudo-cubic. The formation of a lattice superstructure 
may stabilize the magnetic tendency to phase separate and 
minimize lattice distortions. 

Note, however, that phase separation may manifest it- 
self as in "frustrated phase separation" scenarios |fl5"| : 
Since the Coulombic interaction between holes was not 
explicitly included it is possible that in realistic situa- 
tions phase separation may be replaced by the forma- 
tion of complex structures, such as the stripes observed 
in cuprates [^6[ [jj], Thus, it is reasonable to specu- 



late that these stripes could also appear in the insulating 
regime of the manganites and they should be detectable 
using neutron scattering techniques. 

On the theoretical side, future work will be directed to 
the analysis of the influence of phonons and orbital de- 
generacy into the phase diagram observed in the present 
paper, specially regarding phase separation, as well as the 
calculation of dynamical properties for the models investi- 
gated here. Work is in progress along these fronts to com- 
plete a qualitative understanding of the phase diagram 
corresponding to models for the manganites beyond the 
double exchange model. 
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